E.M. Ovsiyuk* 
Particle in the Electromagnetic Wave with 
Cylindrical Symmetry, an Analog of the Volkov 

Problem 

Mozyr State Pedagogical University named after LP. Shamyakin, Belarus 



Abstract 

Cylindrically symmetric analogue of the Volkov problem is exam- 
ined. In presence of the field of a cylindrical electromagnetic wave, 
classical motion of a non-relativistic particle on a cylindrical surface is 
described exactly in terms of elliptic functions. 



1. Particle in a plane wave, the Volkov problem 

Problems of mathematical physics constitute a basis for simulating processes 
in modern physics. Expanding the range of such problems is useful in many 
ways. In particular, in the context of investigating the properties of nanos- 
tructures may be of interest the transport of electrical charges in a very 
narrow cylindrical conductive layer in presence of an external electromag- 
netic wave with cylindrical symmetry. As shown in this study, such a system 
has many common mathematical properties of the well-known Volkov prob- 
lem [1] on the motion of electric charge in the ordinary plane electromagnetic 
wave. 

Let us start with the potential for a plane wave propagating along posi- 
tive direction of the z-axis 

Ai(x^) = Acos—{t- — ) = ^cos(a;t - kx^) , 
T c 
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^0 = 0,^2 = 0,^3 = 0,^=-; 

c 

Foi = -Ak sin(wt - kx^) , F31 = +Ak sm{u)t - kx^) . (1.1) 
Transforming this field to the cyhndrical coordinates 

x^ = r cos (j) , x^ = r sin (j) , x^ = z , 

dS^ = c^dt^ - dr^ - ? d^? - dz^ , (1.2a) 

we get 

Qx\ 

Aj. = -^Ai = +A cos(bcos(ujt — kz) , ^0 = , 
or 

dx^ 

A^ = -^Ai = -Arsm(pcos{Lot-kz) , A^ = ; (1.26) 

Fq^ = ~ Ak cos ([) sm{ujt — kz) , Fq^ = +Ak rcos(/)sin(a;t — kz) , Fq^ = , 
F(t>z = kA r sin ^ sin(wt — kz) , Fzr = kA cos ^ sin(wt — kz) , = . 

(1.2c) 

Next we use the shortening notation for the phase: ujt — kz = Q; then a 
plane wave is described by the relations 

Aj. = +A cos (f) cos O , A^ = — A r sin cos O , 

= — ^A;cos<^ sinQ , Fq^ = +Ak r sin cf) sinU , 

F(f)z = kA rsincf) sinJ7 , F^r = kA cos^ sinfi . (1.3) 

Let us consider the problem of a particle in the plane wave, using La- 
grangian formalism (from the very beginning specifying the problem in cylin- 
drical coordinates) [2] 

L='^ i-gikV'V'') - - 9^kA'V'' = 
= ^ ( v^'V + V'^V't' + V^V^ )- 

--A cos 6 cosO V + -Arsind) cosO . (1.4) 
c c 

Euler-Lagrange equations take the form 

d eA eA 

— ( y cos(/) cosfi ) = rV^ H sin ^ cos 

dt mc mc 



d cA. cA 

— ( r'^V'^ H r sin ^ cos O ) = — sin cos ilV^ H r cos cos O , 

dt mc mc mc 



d cA cA 

—V^ = k cos (I) sinfi H krsincf) sin^^ . 

dt mc mc 



(1.5) 



The system (1.5) is equivalent to 

rV't'V'^ = k cos (t) smn + — kcoscj) sin O , 

dt m mc 

— — + - VV'I' = — k sinJ7 k sinQ , 

dt r m r mc r 

(lyz pj^ 

—— = k coscf) sinnF'' + — krsincf) sin^^y^. (1.6) 

dt mc mc 

A simpler system of equations can be obtained by starting with Cartesian 

coordinates 

dV^ 

m 



m- 



dt 

dV^ 
dt 





- V^B'^) 


= eE^ + e{V^B^ 


- V^B^) 


= eE^ + e{V^B^ 


- V^B^) 



= eE' + e{V'B' - V'B') . (1.7) 
Taking into account 

E'^ = Foi = -Ak sm{ujt - kx^) , 

cB^ = -F31 = -Ak sm{Lot - kx^) , 
equations (1.7) are written 

— — = k smQ-\ k y smCl , 

dt m mc 

—-=0, — - = fcF^sini^. (1.8) 

dt dt mc 

Systems (1.8) and (1.6) must be equivalent, which is readily verified by 
direct recalculation, using the formulas 

= -rsin<^y'^ + cos(^y'' , = r COS + sinewy'' . (1.9) 
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Next we will analyze the equations (1.8). According to (1.8), along the 
•3? the particle moves with constant velocity. In the plane of the 1— 3, its 
motion is described by equations 

; — = K smi2 H k — smSZ , 

c at mc mc c 

— = sinJ^ — . (1.10) 

c at mc c 

Let us introduce the variables 

q= k, [q] = , v^ = — , v^ = — 1.11) 

mc sec c c 

then equations (1.10) can be written as 

—r— = —q sinQ + qv^ sinQ , — — = —q sinO . (1-12) 

You can find an approximate solution of this equation, considering the 
motion for sufficiently small time intervals (see [3]): 

v^«l, kz = - —dt = oj / — dt«ojt (1.13) 
c Jo at J c 

that is, by imposing an additional constraint = ut — kz ^ ut . This system 
of equations (1.12) takes the form 

du^ dt)^ 

—— = —q suiujt , — — = —q sinwi . (1.14) 
dt dt ^ ' 

Integrating the first equation, we get 

V = — cos Lot + {V/r,\ ) , 

- = ^ smut + [vL. )t-\ . (1-15) 

C ^ ' id c 

After that, it is easy to integrate the second equation 



o2 



dv^ r q , , / 1 X 1 . , 

— — = —q — cos cut + (V/r.) j smwr , 

dt oj ^ ' u 

• 2 4. . 1 / 1 1 \ j.ir3 
^ sm ujt + -{v^,^ - -) cosut + [ v^,^ - -(t;(o) 
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z t sm2a;t , Q , i q. . ,,r3 1^1 Q ^ , . 

c = - 2^ ^ 2 - + ^» --)s.nut+[ .(0) - - (^(0) " 1 ^ + ^ • 

(1.16) 

Let us turn to the system (1.12) in the general case, and translate it to 
a new time variable r 

z z dz 
T = t , t = T + - , dt = dT-\ . (1.17a) 

c c c 

Generalized velocities transform according to 

I c~^dx^ 2 c~^dx'^ ""^ 



V 

V" = T-r- = , V 



dr + c ^dz 1 + dr + c ^dz 1 + ^ 

,1 1 dx .n 1 dx , Q 1 dx 

v^ = c ^— , v^ = c ^— , v'^ = c ^— , 
dr dr dr 

c^^dz 

v = -■, TT- = --, ^ ' ^ = --, ? ; (1.176) 

dr + c-^dz l + v^ ^ ' 

operator of differentiation with respect to time is converted according to 

Jt^ " dr + c-'^dz^ " ■ 
The system of equations (1.12) takes the form 

1 d . . 

-, , ^3 TTT^ = ~^ ^^^^'^ + 5 smwT -—To , 
1 + dr 1 + 1 + 

1 d . 



1 + dr 1 + 1 + & 

The previous equation can be rewritten as 

1 Id. -v^ . _ 



—q sinwT ^ . (1-18) 



q'sinwT 1 + dr 1 + l + i)^ 1 + v 

1 1 d , 
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q sin LOT 1 + dr 1 + 1 + 

Allowing for the second equation, from the first one we get 

1 1 d \ / 1 1 d \ 



gsinwr 1 + dr J \ gsinwr 1 + dr J 1 + 1 + v' 



(1.21a) 



or 
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d_ 



1 



1 



d 



sin UT I + dr 1 + 
After integrating over r we have 



q smwr . 



1 



d 



1 + dr 1 + 



( coswT + A) sinwr. 



(1.216) 



(1.21c) 



Note that in the hmit « 1 , kz « ut , t ^ t , ii^ k, , the previous 
equation is simphfied 



d n , q . 
—V = ( cosujt + A) smuit , 

dt ijj 

which coincides with the first equation in (1.16). 
We return to equation (1.21c): 



(1.22) 



1 



d 



1 + t)^ dr 1 + 



( coswT + A) smwT. 



It is a differential equation of first order with separated variables 



d 



{1 + v 



3^2 



( coswT + A) sinwr dr , 



(1.23a) 



(1.236) 



which is easily integrated (we introduce a constant integration) 



1 q^ sm^wT A A 

■ 7rr. r5T^ = 7: cos cut — — „ 

2(l + fi3)2 a; 2a; oj lup- 



and further 



OJ 



\J q^ SIT? ojT + 2Aa; cos wr + A 



(1.23c) 



(1.23d) 



Prom (1.23(i) one can obtain an expression for the particle velocity 
see (1.176): 



1 

c dt 



l + {;3 ' 



dz c I 

— = c ^ — \lq^ sin^ uj{t - zlc) + 2\uj cos uj{t - z/c) + k . (1.24o) 
dt OJ ^ 

If we assume that at the time t = coordinate z = = 0, then from 
(1.24a) it follows 

io = c T - V2Aw + A . (1.246) 

OJ 



Variant with the lower sign must be discarded as not physical, because 
it involves an initial velocity greater than the rate of of light. We return to 
equation (1.23c?) 




— = -1± 



dz 



\J s\v? ujT + 2\u cos wr + A 



so we get 



u = UJT = uj{t — z/c) 




I 



du 



(1.25) 



^Jq'^ + K — cos^ u + 2\uj cos u 



This integral reduces to an elliptic one (in more detail we consider the 
analogous integral in the next section). As a result we get the function of 
z{t) in an implicit transcendental form - this is a known property of the 
classical problem of motion of a charged particle in a plane electromagnetic 
wave (Volkov [1]). 

In the next section we will construct a cylindrically symmetric analog 
of such a system, treating a particle in the simplest plane electromagnetic 
wave with cylindrical symmetry. 

2. Particle in the presence of a cylindrically 
symmetric electromagnetic wave 

We start from a simple solution of the Maxwell's equations in cylindrical 
coordinates {ct,r,(f),z) 

^0 = 0, Ar = 0, A^ = a cos{ujt - kz) , = ; (2.1) 

it meets the following electromagnetic tensor 

Fq^ = —ka sin{ujt — kz) , = +ka sin{ujt — kz) . (2-2) 

It is easy to check that here indeed we have the solution of the Maxwell's 
equations 



1 



d 



^/^F'^^ = 




which is reduced to one non-trivial equation 

^ = ^ , -c5t^Fo0 + 5^^F^<^ = 
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— c ^dt[—ka sm{ajt — kz)] + dz[ka sm{ujt — kz)] = , 
that is an identity 

+k^a cos{uit — kz) — k^a cos{uit — kz) = . 

It is useful to convert the electromagnetic 4-potential to Cartesian coor- 
dinates 

^2 = + (^1)2^(^2)2 « cos(a;i - kz) . (2.3) 
It corresponds to the electromagnetic tensor 

cB' = F23 = - (^1)2 + (^2)2 ^« si^M - kz) , 
cB'^ = F31 = --^-TT^ -r^Ty^ ka sm{cot - kz) , 



(x^)2 + (x 

" " (Xl)2 + (x2)2 ~ ^"""^ ' 

cS^ = F12 = , ^3 ^ ^ . (2.4) 

This is a wave propagating along the axis z, it is singular on the axis z 
(when r = 0). In the plane of the fixed value of the variable z, at each point 
the magnetic field is directed toward the center (0, 0) along the radius, and 
the electric field at each point in the plane is directed along the vector e<^. 
Amplitude of the electric and magnetic fields varies according to the same 
law ^ ^ 

E ^ - ak sin(a;f — kz) , cB ^ - ak sin(a;t — kz) . 
r r 

Consider a particle in that field. The Lagrangian of the system is given 

by (examining the problem of a non-relativistic particle) 

^ = f ( -9r3{^)y'y') - I Qrji^) A'V^ = 

= ^ [v^v'' + r^Vt-V^ + V^V^) - - a cos(a;t - kz)V'(' . (2.5) 
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Euler-Lagrange equations 

d dL dL d dL dL d dL dL 



dt dV^ dr ' dt dV't> d(f) ' dt dV^ dz ' 
here take the form 

dt ' 

[ - — cos(a;t -kz)]=0 =^ ^ = , 
dt^ mc ^ ' ^ dt ' 

^ = -—ksm(ujt-kz)V'^. (2.6) 
dt mc ^ ' ^ ' 

Prom (2.6) we get a much more simpler systems, if consider the motion 

at sufficiently small time intervals: 



CO dz ^ fV^ ^ „ , 

kz = - — dt = uj / — dt «ut (2.7a) 
c Jo dt J c 

that is, by imposing an additional constraint 

uit — kz ^ ujt . (2-76) 
Then (2.6) gives (we use the notation ea/mc = b) 



dt ' ' ' r2 
dV^ 



= -bk sinujtV'f' . (2.8) 

and then we get three equations with separated variables: 

d'^r _ {I + b cos ut)"^ 
dt^ ' 

dV^ _ sincjt (/ + 6 cosoot) 

Unfortunately, such a system hardly can be solved as well. One can try 
to simplify the original system by imposing additional constraint. The most 

interesting from a physical point of view, is the condition r = vq = const. 
It means that a particle can only move along a given cylindrical surface. In 
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this case the electric component of the external wave speed varies linearly 
along the vector e<^, and the magnetic field changes the linear velocity along 
the axis z. 

With this condition, instead of (2.6) we will have a simpler Lagrangian 
and two new equations of motion: 

L = ^ {rlV^yi' + V^V'') - - a cos(a;t - kz)V'^ , 

^ [ rlV^ - b cos{ujt -kz)] = 0, 
fjyz 

= -bk sm{ujt - kz) V'>' . (2.10) 
Variable (j) varies with time according to 



ytf, _ I + b cos{uit — kz) 



(2.11a) 

^0 

The second equation in (2.10) takes the form differential equation for z{t) 

— ; — = sin(a;t — kz) [ I + b cosiut — kz) 1 . (2.116) 

dt Vq 

First wc will find the solution of equations (2.11) in the approximation 
of the velocity much smaller than the speed of light, and not too large time 
intervals 

uj dz , fV^ , 

kz = — — dt = ui / — dt « ujt , ujt — kz uit . 

c Jo dt J c 

In this case (2.11) can be simplified 

<^ _ J + 6 cos cot _ I + b 

= 2 sinoot {I + b coscot) . (2-12) 



dt 

These equations are easily integrated: 



(p{t) = 4>o -\ — 2 ^ "I — 2~ sinwt , (2.13a) 



bk , b . o I . bkl. 



= [ — sin^wt - - coswt 1 + (Vri „ 
ri 2lu uj rico 
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, , f , hk , h 1 — cos 2ujt I , ,^ hkl , , , 

<t) = / {-^ o 7, cosa;t + Fo - —)) dt = 

7 Tq 2a; 2 ui TqUJ 

bk , b ,1 sin2a;t, / . , 

(2.136) 

Character of the motion is as follows: along the z axis the particle moves 
with a certain average constant speed, with two superimposed oscillatory 
movement in time with frequencies lo and 2a;; on variable (p is also the 
movement at a constant (angular) velocity and when imposed oscillating 
motion. 

Let us return to the analysis of equation (2.116) without approximation 
of small intervals of time. It is convenient to introduce a new time variable 

T = t--, t = T + - , dt = dT + — . (2.14a) 
c c c 

Generalized velocity is transformed according to 

yz ^ ^ yz ^ y (2.146) 

dT + c-^dz 1 + c-iy^ 1 - c-iy^ ^ ^ 

operator of differentiation with respect to time is converted according to 

d , _ d 1 d (o-\A) 

dt-^ ~ dr + c-^dz'' ~ l + c-^V^ dr"^ ' ^' ^' 

Equation (2.116) takes the form 

1 d c~^V^ bk sm. (jjT {I + b cos ujt) 

1 + c-iy^ dT 1 + c-iF^ ~ c rl 

or (go to the dimensionless velocity = c~^V^) 

f I \ 2huj 

d 7" TTTT = -^r-n sina;r (/ + 6 cosa;T) dr . (2.15o) 

V (1 + 1)^)27 cVg ^ ^ ^ ' 

In fact, this equation (up to notation) coincides with (1.236) - an equa- 
tion arising in the case of particles in the field of an ordinary plane wave. 
After integration of (2.15a) 

1 26a; . 6 . 2 I 1 ^ 

-, TTT = Sin UT cos UJT \ + C 

(1 + v^f c^rl ^2oj uj ^ 
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and further 

1 , [2bu} , b 



[ sin^ COT- - coswr ] + C = ± ^tF) • (2-156) 



1 + y c^Fq ^ 2a; lo 

The equation for the velocity can be integrated: 
-1 dz .2 1 



There wih arise an imphcitly defined function z{f) with use of ehiptic 
functions. Indeed, consider the integral 

f dw ,^ ^ 

J = —^^^=^^^= ; (2.16a) 
J V a cos^ + 6 cos w -\- c 

introducing a variable 

du 

— cos w = u , aw = , , 

we get 

J = / , . (2.166) 

J ^/{l-u^^){au^^ +bu + c) ' 

It is convenient to introduce yet another variable trough a linear frac- 
tional transformation 

p X + (T 

u = , 

x + 1 

2 _ (1 - p2)a;2 + 2(1 - p a)x + (1 - a^) 
{1 + x 

{ap^ + bp + c)x'^ + [2apa + b{p + a) + 2c]x + (aa^ + ba + c) 



{l + xf 

The required coefficients are determined by conditions 

1 - pa = , 2a pa + 6 (p + (j) + 2c = , 
_ -(a + c) + V(a + c)2-62 
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-(a + c)- V(a + c)2-62 



As a result, the integral J is the following: 

^ , dx 



= 1 / ^ 

A J - p^x^ 



m^x 



A = 



a + pb + p^c 



Making the change of variables 



m 



2 _ p\p^ a + pb + c) 
a + pb + p^c 



p x 



, k = 



m 



(2.176) 



(2.18) 



we reduce the integral to the canonical form of the elliptic integral: 

1 r dU 



J = 



I 



(p2 + m2) ] VI - f/2 Vl - A:2[72 ' 
Therefore, as a result, we produce solution z{i) in an implicit form 



(2.19) 



z{t) 



-CT zb — 



dU 



^JA (p2 + m^) J Vl - ?72 Vl - A;2[/2 



(2.20) 



Author is grateful to V.M. Redkov and Yu.A. Kurochkin for advice and 
for help in writing this paper. 
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